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We show how to capture a single photon of arbitrary temporal shape with one atom coupled to 
an optical cavity. Our model applies to Raman transitions in three-level atoms with one branch of 
the transition controlled by a (classical) laser pulse, and the other coupled to the cavity Photons 
impinging on the cavity normally exhibit partial reflection, transmission, and/or absorption by the 
atom. Only a control pulse of suitable temporal shape ensures impedance matching throughout 
the pulse, which is necessary for complete state mapping from photon to atom. For most possible 
photon shapes, we derive an unambiguous analytic expression for the shape of this control pulse, 
and we discuss how this relates to a quantum memory. 

PACS numbers: 03.67.-a, 32.80.Qk, 42.50.Dv, 42.50.Pq, 42.50.Ex, 42.65.Dr 



In their seminal paper on the quantum internet, 
J. 1. Cirac et al. [1] introduced a scheme for photon- 
mediated state mapping between two distant atoms 
placed in spatially separated optical cavities. The exper- 
imental attainment of this would have important conse- 
quences for the fields of quantum information process- 
ing (QIP) [2, 3], quantum networking, and the ultimate 
realisation of a "quantum internet" [4], where a single 
network link involves the state mapping from atom to 
photon in a first cavity, and then back to another atom 
in a second cavity, as illustrated in fig. 1. So far, good 
progress has been made experimentally: entanglement 
mapping mediated by photons with atomic ensembles 
[5, 6], cavity QED atom-photon state mapping [7-9], 
single-atom single-photon absorption in ion traps [10] and 
single-photon absorption by a coupled atom-cavity sys- 
tem [11, 12] have all been demonstrated. However, no 
one has yet managed to combine all of the elements nec- 
essary to implement the scheme in its entirety. Further, 
there exists no general method for determining the opti- 
mal control of the classical laser pulses used in the gener- 
ation and absorption of single photons in these processes, 
although some work has been done towards this [13-18]. 
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FIG. 1. (Color online) The Quantum Internet, showing the 
state of one atom (stationary qubit) being transferred via a 
photon (flying qubit) to that of another atom in a separate 
cavity. 

Based on our recent discussion and demonstration of 
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FIG. 2. (Color online) Coupling schemes. The left schematic 
shows the input-output cavity coupling formulation, whilst 
the right diagram shows the atomic states |e), \x) and \g) 
with their respective cavity photon number states |0) and |1). 
The couplings of the control pulse and cavity are given by 
Q(t) and g respectively. 



the single-photon generation process in the time domain 
[19, 20], we here devise an analytical method for finding 
the control pulse required for the complete absorption of 
single photons of arbitrary temporal shape, i.e. for any 
given running- wave probability amplitude (j)in{t) imping- 
ing on a cavity mirror [21]. We emphasise that this tech- 
nique applies directly to the implementation of a quan- 
tum memory and is pertinent to a variety of cavity-based 
systems. As will be discussed later, this scheme relates 
most obviously to mapping Fock-state encoded qubits to 
atomic states [1, 9], but also very simply extends to other 
possible superposition states, e.g. photonic time bin or 
polarisation encoded qubits [7, 12, 22]. 

In the original transmission protocol [1], the leading 
and trailing edges of the photon change role upon its ab- 
sorption by the second atom. The claim had therefore 
been made that such an approach would require photons 
of symmetric temporal shape; it has become evident that 
this is not a strict condition, and models have been made 
using a frequency-mode decomposition of free space to il- 
lustrate this mathematically [15, 23]. Our procedure op- 
erates in the time domain and thus does not require any 
such decomposition into frequency modes. It is there- 
fore closely analogous to the photon generation process. 
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where the temporal shape of an emitted photon is unam- 
biguously linked to the shape of the driving laser pulse. 

Prior to investigating the effect of the atom-cavity and 
atom-laser coupling, we shall briefly formulate a descrip- 
tion of the input-output coupling of an optical cavity in 
the time domain. Inside the cavity, we assume that the 
mode spacing is so large that only one single frequency 
mode is involved, with the dimension-less probability am- 
plitude Ccav{t) for occupying the one-photon Fock state 
|1) at resonance frequency uJcav Furthermore, we assume 
that one cavity mirror has a reflectivity of 100 %, thus 
ensuring that coupling to the outside field is controlled 
uniquely by the field refiection and transmission coeffi- 
cients, r and r, of the other mirror. It is then very con- 
venient to decompose the cavity mode into submodes |-h) 
and |— ), travelling towards and away from this mirror, 
respectively. In turn, the spatio-temporal representation 
of the cavity field reads 



4>^(t)\+) + 4>-{t)\-), 



(1) 



with = + A^, where any change of the 

running-wave probability amplitude from ^+ to 0_ at 
the mirror is taken into account with Acf). The latter is 
small for mirrors of high refiectivity, such that Ccav (t) — 
(j)j^{t)^/t^ :^ (j)-{t)^/t^^ and = 2L/c the cavity's round- 
trip time (note that all running-wave probability ampli- 
tudes (j)any ctrc of dimension s~^^'^ ^ whereas the occupa- 
tional probability amplitudes Cany are dimensionless) . To 
properly describe the field outside the cavity, we decom- 
pose it also into incoming and outgoing spatio-temporal 
field modes, with running- wave probability amplitudes 
0m + and (j)out{t — zjc) for finding the photon in 
the I in) and \out) states at time t and position 2:, respec- 
tively. The coupling mirror at z = acts as a beam split- 
ter with the operator aL(ra+ -h rain) + aJ^^(Ta+ — rain) 
coupling the four running- wave modes inside and outside 
the cavity. In matrix form, this coupling equation reads 



r r 
T — r 



(2) 



To relate Ccavify to the running- wave probability ampli- 
tudes, we take r 1 — kI^ and r = \/2/^t^, where k is 
the polarisation decay rate of the cavity. Furthermore, 
we make use of 

^C,cav ^ Ccav{t-\-tr) Ccav{t) ^ ((/)_ 0+)v^ ^ 

dt tip tip trp 

(3) 

With these relations, the first line of Eq. (2) can be writ- 
ten as Ccdy/tj- ~\~ Ccdy ~ (1 f^tr^c-Qdy /tj- H" Tficrc- 
fore Eq. (2) takes the form of a differential equation 
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which describes the coupling of a resonant photon into 
and out of the cavity mode. We emphasise that any de- 



viation from resonance with the cavity mode can be in- 
cluded in time-dependent phase factors of the probability 
amplitudes. Therefore the above coupling model applies 
to any case where just one cavity-field mode is involved, 
using the frequency Uc of this mode as a carrier. We 
would like to emphasise that the result obtained from our 
simplified input-output model is fully equivalent to the 
conclusions drawn from the more sophisticated standard 
approach that involve a decomposition of the continuum 
into a large number of frequency modes [24]. No such 
decomposition is required here as we model the system 
uniquely in the time domain. Note also that the expan- 
sion of the present model to both mirrors being transpar- 
ent is straight forward, but not required for the following 
analysis. 

Next, we need to examine the coupling of a single atom 
to the field mode of the cavity. This has been discussed in 
a recent review [25] and the references mentioned therein. 
To best illustrate the process, we consider a three level A- 
type atom with two different electronically stable ground 
states |e) and \g), coupled by either the cavity field mode 
or the control laser field to one-and-the-same electron- 
ically excited state \x). For the one-photon multiplet 
of the generalised Jaynes-Cummings ladder, the cavity- 
mediated coupling between 1) and |x, 0) is given by the 
atom-cavity coupling strength while the control laser 
couples |e,0) with |x,0) with Rabi frequency The 
probability amplitudes of these particular three product 
states read Ce(t), Cg{t), and Cx{t), respectively, with their 
time evolution given by the master equation of the cou- 
pled atom-cavity system 



'out ) 



( 



\ 










-iO(t)72 
-7 







-K 







2Ki 







) 


\(t>in J 



(5) 

which combines the Schrodinger equation for a three-level 
atom coupled to the cavity with the input-output relation 
from Eq. (4). Here 7 is the spontaneous-emission rate of 
the atom in state x. The coupling to the external field is 
taken into account by the mirror- induced decay rate k,. 

Note that Cg{t) = Ccav{t) because the state |^, 1) is the 
only atom-field product state in which there is one pho- 
ton in the cavity. We now restrict our discussion to the 
resonant case with the overall phase chosen such that g 
and (j)in{t) are both real for all t. By consequence, the 
probability amplitudes Cg{t) and Ce{t) are real as well, 
while Cx{t) is purely imaginary. We also emphasize that 
we are considering only one quantum of excitation, there- 
fore at any given time, the probability of occupying 0) 
is given by the total probability of having a photon out- 
side the cavity, i.e. either in state \in) or in \out). The 
only way of coupling these states is via the cavity mirror 
to 1^, 1), which is taken into account by Eq. (5). 

Our goal is to completely absorb an incoming photon. 
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with its wavepacket given by <j)in{t). Obviously, com- 
plete impedance matching is a necessary condition, i.e. 
(l>out{t) = must be met at all times. We furthermore 
assume that r ~ 1, i.e. the reflectivity of the cavity mir- 
ror is nearly one. Using these constraints with eqn. (5) 
we express two of the probability amplitudes in terms of 
the specified photon function ^in(^) and its time deriva- 
tive, 



(6) 



Cg{t) = (l)in{t)/V2K 



= I 



(t>in{t) - i^(t>in{t) /g*y2n. 



(7) 



We also obtain the product 



^{t)ce{t) = 2 [zc,(t) + Z7c,(t) - gcgit)] = C{t). (8) 

To proceed, we need to consider the initial population 
of states and the overall continuity of probabilities (by 
which we mean the excitation is conserved). It is nat- 
ural to assume that the atom-cavity system is initially 
prepared in state \g, 0), which lies outside the considered 
subspace. The photon is completely in the incoming state 
|m), i.e. / \ (pin{t)\'^dt = 1. The portion that couples into 
the cavity therefore directly populates 1), so that the 
continuity balance yields 

Pee{t) = P0-Pgg{t)-P:rx{t)^ [ [\cl)in{t')\^ -^2^ Pxx{t')W . 

(9) 

Here, pa = 0*01 are the populations and therefore the 
diagonal elements of the density matrix, and 2^pxx is 
the population decay rate from the electronically excited 
state. To account for an imperfect state preparation with 
a small initial population in state |e, 0), the offset term po 
has been introduced phenomenologically. The relevance 
of this term (which is ideally zero) becomes obvious in 
the following discussion. 

Equation (9) now gives a direct analytical expression 
for Pee 7 and as Cg is real on resonance, we simply divide 
Eq. (8) by \/~p^J^ to obtain 



C{t) _ 2[ic^{t)^i-fc^{t)-gcg{t)] 



\/peeit) 



(10) 



This is an analytical expression, derived following equa- 
tions (6) to (9), for the Rabi frequency ^}{t) required to 
achieve full impedance matching over all times and there- 
fore to absorb the incoming photonic wave packet (pinit) 
completely by the atom-cavity system. Most computer 
algebra systems (CAS) can be used to obtain a closed ex- 
pression for rt{t) given the a functional expression for the 
incoming photon. Of course, this non-iterative algorithm 
may also be applied numerically. 

We now consider physically realistic photons that are 
normally restricted to a finite support of well-defined 
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FIG. 3. (Color online) Incoming sin^ photon (dotted). Case 
(a): Empty cavity, all reflected (dashed); Case (b): System 
prepared in |^,0), small reflection (dash-dotted); Case (c): 
Small initial population in |e, 0), reflection suppressed (thin 
solid) . The control pulse (thick solid) is derived to match case 
(c). 
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FIG. 4. (Color online) A twin peak pulse. All traces as de- 
scribed in fig. 3. 



start and end times, t start and tstop- We also assume they 

start off smoothly, i.e. with (l)in{tstart) = ^ 4> in {t start) = 

0, as described in [19]. However, the second derivative 
might be non-zero at t starts so that Eqn. (8) yields 
^{t start) Ce{t start) 7^ 0. To Satisfy the latter inequality 
a small but non-vanishing initial population is required 
in the state |e,0). In other words, perfect impedance 
matching with po = would only be possible with pho- 
tons of a physically impossible infinite duration. 

To illustrate the power of the procedure and the im- 
plications of the constraints to the initial population, 
we now apply the scheme to a couple of typical photon 
shapes that one may obtain from atom-cavity systems. 
To this purpose, we consider a cavity with parameters 
similar to one of our own experimental implementations. 
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with 7) = 27r X (15, 3, 3) MHz and a resonator length 
of L = 100 fim. As a first example, we assume that a sym- 
metric photon with (j)in{t) cx sin^ {nt / Tphoton) impinges on 
the cavity. For a photon duration of Tphoton = 3.14 /is, 
fig. 3 shows (j)in{i)^ ^{t) and the probability amplitude of 
the reflected photon, (poutif)^ as a function of time. The 
latter is obtained from a numerical solution of eqn. (5) 
for the three cases of (a) an empty cavity, (b) an atom 
coupled to the cavity initially prepared in |^,0), with 
Po = 0, and (c) a small fraction of the atomic population 
initially in state |e, 0), with po = 0.5%. In all three cases, 
the Rabi frequency Vt(t) of the control pulse is identical. 
It has been derived analytically assuming a small value 
of po = 0.5% (this choice is arbitrary and only limited by 
practical considerations, as will be discussed later). From 
these simulations, it is obvious the photon gets fully re- 
flected if no atom is present (case a), albeit with a slight 
retardation due to the finite cavity build-up time. Be- 
cause the direct reflection of the coupling mirror is in 
phase with the incoming photon and the light from the 
cavity coupled through that mirror is out-of phase by tt, 
the phase of the reflected photon flips around as soon as 
Ccav{t) = <l^in{t)/V^' This shows up in the logarithmic 
plot as a sharp kink in (j)out(f) around t = 0.13 /is. 

The situation changes dramatically if there is an atom 
coupled to the cavity mode. For instance, with the ini- 
tial population matching the starting conditions used to 
derive Vt{t)^ i.e. case (c) with po = 0.5%, no photon 
is reflected. The amplitude of \(j)out{t)\^ remains below 
10~^^, which corresponds to zero within the numerical 
precision. However, for the more realistic case (b) of the 
atom-cavity system well prepared in |^, 0), the same con- 
trol pulse is not as efficient, and the photon is reflected 
off the cavity with an overall probability of 0.5%. This 
matches the "defect" in the initial state preparation, and 
can be explained by the finite cavity build-up time lead- 
ing to an impedance mismatch in the onset of the pulse. 

We emphasise that this seemingly small deficiency in 
the photon absorption might become significant with 
photons of much shorter duration. For instance, in the 
extreme case of a photon duration Tphoton < build- 
ing up the field in the cavity to counterbalance the direct 
reflection by means of destructive interference is achieved 
most rapidly without any atom. Any atom in the cavity 
will act as a sink, removing intra-cavity photons. With 
an atom present, a possible alternative is to start off with 
a very strong initial Rabi frequency of the control pulse. 
This will project the atom-cavity system initially into a 
dark state, so that the atom does not deplete the cavity 
mode. Nonetheless, the initial reflection losses would still 
be as high as for an empty cavity. 

Second, we consider a more sophisticated non- 
symmetric twin-peaked photon impinging onto the 
cavity, with (l)in{t) cx sin^(27rt/rp/,oton) cos(7r/2(l - 
t / Tphoton)) • As in the first example, we simulate the pro- 
cess with a control-pulse Rabi frequency Q{t) derived for 




FIG. 5. (Color online) Derived Rabi frequency, ^^(t), as a 
function of t and po. In this particular case, the required 
pulse-shape changes little for values of po greater than 0.002. 

po = 0.5%. From fig. 4, it is obvious the results are the 
same as before, despite the rather complicated photon 
shape. This impressively demonstrates the potential of 
our method, giving one a very simple procedure to de- 
rive control-laser pulses for absorbing incoming photons 
of arbitrary shape. 

Finally, we investigate how Q{t) changes in the limit 
po — > 0, which would normally correspond to the ini- 
tial condition in any implementation. For the twin-peak 
photons discussed above, fig. 5 shows Q{t) for different 
values of po- As one would expect from our above dis- 
cussion of the build-up dynamics, Q{t) does not converge 
to a finite function in the limit po — > 0. Hence for any 
realistic scenario, one needs to chose po to be as small 
as possible such that a feasible driving pulse shape is ob- 
tained. Limits to this are normally the finite laser power 
and/or the finite bandwidth of the amplitude modula- 
tors. For the 'best possible' driving pulse derived in this 
manner, a numerical simulation starting with the atom- 
cavity system prepared in |^, 0) then reveals the actual 
efficiency. In the particular examples discussed here, the 
photon-reflection losses are comparable to the non-zero 
Po used to derive the driving pulse, i.e. at a level below 
0.5% which is negligible compared to the noise affecting 
current experimental approaches. 

A unique feature of the method presented here is its 
universal applicability to a wide range of atom-cavity 
coupling regimes. The simple recipe for calculating the 
control pulse assures full impedance matching for any 
given single-photon wavepacket impinging on the cav- 
ity. The present work therefore extends the recent ap- 
proach made by Gorshkov et al. [17], as we need no a 
priori restriction to extreme coupling cases in their re- 
spective approximations. While the authors of the for- 
mer were seeking for the best possible quantum memory. 
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FIG. 7. (Color online) Mapping of a photonic time-bin super- 
position state to a superposition of atomic (spin) states. The 
left diagram shows a photonic qubit impinging on an atom- 
cavity system, and the respective driving pulses required in 
order to map the qubit to the atomic states (right). 



FIG. 6. Storage efficiency (probability of mapping the in- 
coming photon to |e, 0), solid line) and impedance mismatch 
(probability of back-reflecting the incoming photon, dashed 
line) as a function of the cooperativity, C = g'^/{2Hij). The 
dotted line shows the maximum possible efficiency. Both ef- 
ficiency and impedance mismatch have been numerically cal- 
culated for At = 7 = 27r3 MHz using symmetric sin^ pulses of 
3.14 /is duration, with control pulses modelled according to 
Eqs.(7-10) to achieve best impedance matching. 



we here ask for perfect impedance matching. The lat- 
ter is a necessary condition (satisfied by one and only 
one control pulse), which may not ensure 100% efficiency 
under all circumstances. To illustrate the interplay of 
impedance matching and memory efficiency, Fig. 6 shows 
the reflection probability and the memory efficiency (ex- 
citation transfer to |e, 0)) as a function of the cooperativ- 
ity, C = /{2i<i^). Obviously, the impedance matching 
condition is always met, but the efficiency varies. For 
C > 1, it asymptotically reaches the predicted optimum 
[17] of 2C/(2C + 1), but it drops to zero at C = 1/2 (i.e. 
ioT g = n = '^). In this particular case, the spontaneous 
emission loss via the atom equals the transmission of the 
coupling mirror. Hence the coupled atom-cavity system 
behaves like a balanced Fabry- Perot cavity, with one real 
mirror being the input coupler, and the spontaneously 
emitting atom acting as output coupler. Therefore the 
photon gets not mapped to |e,0). This limiting case fur- 
thermore implies that impedance matching is not pos- 
sible for C < 1/2, as the spontaneous emission via the 
atom would then outweigh the transmission of the cou- 
pling mirror. The application of our formalism therefore 
fails in this weak coupling regime (actually, the evalua- 
tion of Eq. (9) would then yield values of pee < 0, which 
is not possible). 

We now discuss the application of our scheme to map- 
ping an arbitrary time-bin superposition photonic state 
[26], (j)ph{t) = a(j)i{t)-\-/3(j)2{t) to a superposition of atomic 
spins, a\m = —1) + /3|m = +1). Assuming a tt polarised 
incoming photon and and (j)2 not overlapping in time, 
(i.e. they represent two orthogonal field modes), we de- 
vise two driving pulses fli{t) and ft2{t) that map these 



two field modes into \m = —1) and \m = +1), respec- 
tively. As illustrated in fig. 7, the latter can be achieved 
with fti{t) and 1^2 (^) having different circular polarisa- 
tions. As a simple consequence of the unitarity of the 
time evolution, successive application of the two driv- 
ing pulses, each within its respective time bin, results 
in the desired mapping of photonic to atomic superposi- 
tion states. We emphasise that the pair of driving pulses 
only depends on the mode functions ^1,2 describing the 
two time bins, and not on their relative amplitudes a 
and p. Further, the three- level systems addressed dur- 
ing the first and second time bin are different from one 
another, and no cross talk occurs due to the different 
polarisations of the driving pulses. Hence the numerical 
simulation depicted in fig. 3 now applies to the two time 
bins individually, with a specific example of the scheme 
shown in fig. 8. The plot is calculated using the same 
values of (^,^^,7) as before, and demonstrates a storage 
efficiency of 95.3% with a fidelity of unity. The latter is a 
consequence of the fact that the only incoherent process 
in the model, non-adiabatic excitation and decay to and 
from the excited state |x,0), is considered as a loss from 
the system. Thus it reduces the efficiency of the map- 
ping, but does not contribute to any incoherence in the 
final atomic state superposition. 

One could also map a superposition of photonic num- 
ber states to a pair of atomic states [1, 9]. For instance, 
a\0) + /3|1) gets trivially mapped to a\g) + /3|e); if no 
photon arrives then the atom remains unaffected (as the 
atom-cavity system stays in 0)), whilst the signed driv- 
ing pulse assures that a photon arriving in mode (j)ph{t) 
gets mapped to |e), in the way we depicted in fig. 3. 
Again, the unitarity of the time evolution assures the 
accurate mapping of superposition states. 

The photon reabsorption scheme discussed here, to- 
gether with the earlier introduced method for generat- 
ing tailored photons [19, 20], constitute the key to ana- 
lytically calculate the optimal driving pulses needed to 
produce and absorb arbitrarily shaped single-photons (of 
finite support) with three level A-type atoms in optical 
cavities. This is a sine qua non condition for the success- 
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FIG. 8. (Color online) Simulation of the photon to atom state- 
mapping utilising the procedure outlined above. The photonic 
qubit is prepared in the state ^(^i — ^2) and mapped onto 
the atomic spin states \m = —1) and \m = +1). The popu- 
lations of the atomic states are shown (thick line: \m = — 1), 
thick dotted line: \m = +1)), as is the incoming photonic 
state (thin dashed lines). The two driving pulses are shown 
below (arb. units), whilst the upper inset shows the atomic 
state density matrix after the successful absorption of the 
photon. 

ful implementation of a quantum network. It is expected 
that this simple analytical method will have significant 
relevance for those striving to achieve atom-photon state 
transfer in C-QED experiments, where low losses and 
high fidelities are of paramount importance. 
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